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Abstract 



m 

' We study the strong and electromagnetic decay properties of scalar mesons above 

1 GeV within a chiral approach. The scalar-isoscalar states are treated as mixed 
states of quarkonia and glueball configurations. A fit to the experimental decay 
rates listed by the Particle Data group is performed to extract phenomenological 
constraints on the nature of the scalar resonances. A comparison to other theoretical 
^ . approaches in the scalar meson sector is discussed. 
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The interpretation of scalar mesons constitutes an unsolved problem of hadro- 
nic QCD. Below the mass scale of 2 GeV various scalar states [1] are encoun- 
tered: the isoscalar resonances a = / (400 - 1200), / (980), / (1370), / (1500) 
and /o(1710), the isovectors ao(980) and cto(1450) and the isodoublets Kq(800) 
and Kg (1430). The existence of the K^(800) is still not well established and 
omitted from the summary tables of [1] . From a theoretical point of view one 
expects the scalar quark- ant iquark ground state nonet ++ , a scalar-isoscalar 
glueball, which lattice QCD predicts to be the lightest gluonic meson with a 
mass between 1.4-1.8 GeV [2], and possibly other exotic states (non qq states), 
like e.g. four quark states or mesonic molecules [3]. Various interpretations of 
and assignments for the physical scalar resonances in terms of the expected 
theoretical states have been proposed (see, for instance, the review papers [3, 4] 
and References therein). 

In this work we follow the original assignment of Ref . [5] , where in a minimal 
scenario the bare quarkonium states N = ^Jl/2(uu+dd) = nn, S = ss and the 
bare scalar glueball G mix, resulting in the three scalar-isoscalar resonances 
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/o(1370), /o(1500) and /o(1710). Such a mixing scheme has been previously 
investigated by many authors, as for example in the lattice study of [6] or 
within the model approaches of [7-11]. The mesons a (1450) and Kq(1A30) 
are considered as the 1 = 1 and 1 = 1/2 quarkonia J PC = ++ states. On the 
other hand, the scalars below 1 GeV could be interpreted as four quark states 
or mesonic molecules. 

In this Letter, starting from an effective chiral Lagrangian, we perform a tree- 
level analysis of the strong and electromagnetic decays of scalar mesons settled 
in the energy range between 1 and 2 GeV. Although a chiral approach cannot 
be rigorously justified at this energy scale, since loop corrections can be large, 
we intend to use this framework as a phenomenological tool to extract possible 
glueball-quarkonia mixing scenarios from the observed decays. This study is 
in line with the original application of Ref. [14], where detailed arguments can 
be found to justify this approach. 

Chiral Lagrangian - The strong and electromagnetic decays of scalar mesons 
are based on the effective chiral Lagrangian C e s, derived in Chiral Perturbation 
Theory (ChPT) [12-14]. The Lagrangian involves the nonets of pseudoscalar 
(V = £§ =0 PX/Vl) and of scalar mesons (<S = ELo^Vv 7 ^), the electro- 
magnetic field and, in addition, a new degree of freedom, the bare glueball 
field G. We display the lowest order Lagrangian C e s , where the starting point 
is the large N c limit. Corrections to the large N c limit are encoded in £^ ix and 
C s - 



: eff = ^ (D,U D»rf + X + ) + \ (D»SD»S - M 2 S S 2 ) + \{d^G - MqG 2 ) 
+ c s d (Su fl un+c s m (S X+ ) +^G(u,u») + ^=G( X+ ) 



c 



■9 



+ < {S F+ F + ^) + ^=G{ F; v F + ^) + £- ix + C ix . (1) 



Here the symbol (...) denotes the trace over flavor matrices. The constants c s d , 
c s m1 c g d , c 9 m1 Cg and cf define the coupling of scalar fields and the bare glueball 
to pseudoscalar mesons and photons, respectively (the coupling constant cf is 
supposed to be suppressed because gluons do not couple directly to photons. 
However, an intermediate state of two vector mesons can in the framework 
of vector meson dominance generate a coupling of the glueball to the two- 
photon final state; this coupling is however not considered in the numerical 
analysis). The field G is introduced as an extra flavor-singlet field in the ef- 
fective lagrangian. We use the standard notation for the basic blocks of the 
ChPT Lagrangian [12]: U = u 2 = exp(iV\^2/ F) is the chiral field collecting 
pseudoscalar fields in the exponential parametrization, denotes the chiral 
and gauge-invariant derivative, w M = iu^D^Uu^, x± — u^x^ ± u X^ u i X — 
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2B (s + ip), s — M. + . . . and F+ u = ( u^F^Qu + uF^Qu^ ) , where F^ u is 
the stress tensor of the electromagnetic field; Q = diag{2/3, —1/3, —1/3} and 
Ai = diag{m, rh, m s } are the charge and the mass matrix of current quarks, 
respectively (we restrict to the isospin symmetry limit with m u = m d = rh); 
B is the quark vacuum condensate parameter. Then the masses of the pseu- 
doscalar mesons in the leading order of the chiral expansion are given by 
Ml = 2mB , = (m + m s )B , Mfr = (2/3) (rh + 2m s )B . 

We intend to employ the Lagrangian (1) for the tree-level calculation of the 
two-pseudoscalar decays of scalar mesons. At the energy scale of interest, 
E ~ Ms ~ 1.5 GeV, a calculation of loops and an application of the power 
counting rules are not rigorously justified. The aim of the present approach 
is a phenomenological study of scalar meson physics, for which a tree-level 
calculation represents a useful analysis of the issue. 

The term C^ ix ~ rfrf gives rise to the singlet-octet mixing in the pseudoscalar 
and the scalar multiplets. The physical states rj and rj' are given by 

7]° = r\ cos Op — rj sin 9 P , if = r\ sin 6 P + rj cos 9 P , 

where dp is the pseudoscalar mixing angle. We follow the standard proce- 
dure [13-15] and diagonalize the corresponding rf-rf mass matrix to obtain 
the masses of r] and rf. By using M n = 139.57 MeV, M K = 493.677 MeV (the 
physical charged pion and kaon masses), M n = 547.75 MeV and M v > = 957.78 
MeV the mixing angle is determined as 9 P = —9.95°, which corresponds to 
the tree-level result (see details in Ref. [15]). Higher order corrections in ChPT 
cause a doubling of the absolute value of the pseudoscalar mixing angle [15]); 
in our work we restrict to the tree-level evaluation, we therefore consistently 
use the corresponding tree-level result of 9 P = —9.95° (in the present ap- 
proach we do not include the neutral pion when considering mixing in the 
pseudoscalar sector, because we work in the isospin limit; this mixing is small, 
and can be safely neglected when studying the decay of scalar resonances into 
two pseudoscalars) . Similarly, for all pseudoscalar mesons we use the unified 
leptonic decay constant F, which is identified with the pion decay constant 
F = F w = 92.4 MeV. A more accurate analysis including higher orders should 
use the individual couplings of the pseudoscalar mesons (for a detailed discus- 
sion see Refs. [16]). 

Glueball- Quarkonia Mixing - We briefly discuss the mixing of the glueball with 
the scalar-isoscalar quarkonia states. The singlet (S°) and octet (S 8 ) scalar 
quarkonia states are defined in terms of the nonstrange iV and strange S 
components with the respective flavor content [uu + dd\/y/2 and ss as 

S° = + y/l/ZS, S 8 = y/l/3N- ^2/35. (2) 
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In (1) the scalar nonet states have the same mass Ms, corresponding to the 
flavor and large N c limits. Deviations from this configuration (i.e. from the 
large N c limit) are encoded in £^ lix including quarkonia-glueball mixing lead- 
ing to different masses for the scalar mesons (in [14] this is explicitly shown 
for the scalar nonet, but no glueball is included; the glueball mixing is also 
a deviation from large N c , since in this limit the glueball and the quarkonia 
states decouple). As a result, the scalar-isoscalar sector can be described by 
the most general Klein-Gordon (KG) Lagrangian including mixing among the 
N, S and the bare glueball G configurations [5-9]: 

£kg = -\ E *P + Ml\Q-fGS-y/2frGN-eNS, (3) 

1 <i>=N,G,S 

where □ = <9 M <9 M . The parameter / is the quarkonia-glueball mixing strength, 
analogous to the parameter z of Refs. [5-9]; z refers to the quantum mechanical 
case where the mass matrix is linear in the bare mass terms, / in turn is 
related to the the quadratic Klein-Gordon case. The connection between / 
and z discussed in Refs. [8, 11] leads to the approximate relation / ~ 2zMq ■ 
If r = 1, then the glueball is flavor blind and mixes only with the quarkonium 
flavor singlet S* (this is the case in Refs. [5,8,9]); a value r ^ 1 takes into 
account a possible deviation from this limit. However, as deduced from the 
analyses of Refs. [6,7, 11] r is believed to be close to unity. In the following 
we will use the limit r = 1, i.e. we restrict to a flavor blind mixing. 

The parameter e induces a direct mixing between N and S quarkonia states. 
This effect is neglected in [5, 7-9], where flavor mixing is considered a higher- 
order effect. However, a substantial N-S mixing in the scalar sector is the 
starting point of the analysis of Refs. [17, 18]. The origin of quarkonia flavor 
mixing is, according to [17, 18], connected to instantons as in the pseudoscalar 
channel, but with opposite sign (see also [19]): the mixed physical fields are 
predicted to be a higher lying state of flavor structure [Ny/2 — S]/\/Z and a 
lower one with [N + Sy/2]/y/3. Here we study the case e ^ 0, more precisely 
e > 0, which leads to the same phase structure as in Ref. [17,18], but the 
quantitative results and interpretation will differ. The physical scalar states 
\i) are identified with i = fi = /o(1370), f 2 = /o(1500) and f 3 = / (1710), 
which are set up as linear combinations of the bare states: \i) = B' tN \N) + 
B iG \G) + B iS \S). The amplitudes are the elements of a matrix B which 
diagonalizes the mass matrix of bare states including mixing, which in turn 
gives rise to the mass matrix of physical states Q! = diag{M| 1 , M| 2 , M| 3 } . 

Results for suppressed glueball decays - In the following we determine a best fit 
of the parameters entering in Eqs. (1) and (3) to the experimental averages of 
masses and decay modes listed in Ref. [1]. Decay rates resulting from Eq. (1) 
are evaluated on the tree-level. First we analyze the case of a non-decaying 
glueball, i.e. c 9 d = c 9 m = cf = 0. In this scheme the decays are dominated 
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by the quarkonia components [5] in line with large N c arguments. The phe- 
nomenological analysis of Ref. [9] confirmed this trend, but the recent fit of [7] 
shows a strong contribution by the direct decays of the glueball configuration. 
The strength of the glueball decays remains a point to clarify. Our results for 
suppressed glueball decays are given in Table 1 in comparison to the data of 
Ref. [1]. The only accepted average not included in the fit is Ff 2 ^ m r for the rea- 
son that the decay channel f]f]' is produced at threshold, therefore a significant 
distortion due to the finite width of the state is expected. The state / (1710) 
has only been observed in the decays into two pseudoscalar mesons [1]. The 
decay into the final state 4n, which can be fed by higher meson resonances, is 
suppressed [20] . We therefore impose the additional condition that the sum of 
partial decay widths into two pseudoscalar mesons (Tf 3 ) 2P saturates the total 
width (Tf a ) tot with (Tf 3 ) 2p = (Tf 3 ) tot as indicated in Table 1. Such a constraint 
is necessary to obtain meaningful total decay widths: without this condition 
on the full width a minimum for x 2 is obtained where (Tf 3 ) tot is larger than 
1 GeV, a clearly unacceptable solution. A local minimum for Xtot is found for 
the set of parameters 

M N = 1.455 GeV, M G = 1.490 GeV, M s = 1.697 GeV, / = 0.065 GeV 2 , 
e = 0.211 GeV 2 , c s d = 8.48 MeV, c s m = 2.59 MeV xL = 29 - 01 • ( 4 ) 

Table 1. Mass and decay properties of scalar f mesons. 



Quantity 


Exp 


Theory 




M h (MeV) 


1350 ±150 


1417 


0.202 


M h (MeV) 


1507 ±5 


1507 


~ 


M h (MeV) 


1714 ±5 


1714 


0.003 


F f2 ^ (MeV) 


38.0 ±4.6 


38.52 


0.011 


r^gK (MeV) 


9.4 ± 1.7 


10.36 


0.322 


Y h ^ m (MeV) 


5.6 ± 1.4 


1.90 


8.109 


r/ 3 — irn/Tf^KK 


0.20 ±0.06 


0.212 


0.036 




0.48 ±0.15 


0.249 


2.446 


1" ao^KK /r \ Q ^Trri 


0.88 ±0.23 


0.838 


0.032 




0.35 ±0.16 


0.288 


0.150 


T K ^ Kn (MeV) 


273 ± 51 


59.10 


17.590 


(r /3 ) 2p (MeV) 


140 ± 10 


143.27 


0.110 


Xtoi 






29.01 
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Implications of the fit - The bare non-strange quarkonia field N has a mass of 
M^r = 1.455 GeV which is, as desired, similar to the scale set by the isotriplet 
combination a (1450) with a mass of M ao = 1.474 ±0.019 GeV [1]. The mass of 
the bare glueball M G = 1.490 GeV is in agreement with the lattice results [2] 
and with the phenomenological analyses of [5,7,9]. The bare state S has a 
mass of M s = 1.697 GeV, which is about ~ 200 MeV heavier than the N 
state, an acceptable mass difference like in the tensor meson nonet. 

For the glueball-quarkonia mixing parameter we get / = 0.065 GeV 2 , which by 
the approximate relation / ~ 2zM G [8, 11] corresponds to z ~ 21.8 MeV. The 
results of Refs. [7, 9, 11] are z = 85 ± 10 MeV, z = 80 MeV and z ~ 62 MeV, 
respectively, i.e. of the same order, but larger. The introduction of additional 
flavor mixing between the quarkonia configurations in the fit, as done here, 
leads to a reduction of the strength parameter /. The lattice result of Ref. [6] 
with 43 ±31 MeV is in agreement with the present evaluation, but has a large 
uncertainty. 

The flavor mixing parameter resulting from the fit is e = 0.211 GeV 2 . In 
the limit / = the mixed physical states are = 0.97 \N) + 0.26 \S) and 
|/ 3 ) = -0.26 |iV) + 0.97 \S) (and, of course, |/ 2 ) = \G)). The phase structure of 
the mixed states is, as discussed previously, as in [17-19]. But here the strength 
of flavor mixing is smaller, resulting in mixed states which are dominantly iV 
or S. The influence however of (an even small) flavor mixing in strong and 
electromagnetic decays may be non-negligible. 

The mixing matrix B relating the physical to the bare states in the present 
fit is expressed as: 



/ j/ 1 > = j/ (1370)>^ 



I/2) = |/o(1500)> 
l/a) = |/o(1710)> 



^0.86 0.45 0.24 ^ 



-0.45 0.89 -0.06 
-0.24 -0.06 0.97 



\N) = \nn) 

\G) = \99) 
\\S) = \ss) J 



(5) 



The physical resonances are dominated by the diagonal bare components, 
qualitatively in line with Refs. [3, 5, 7, 9]. Since the glueball does not contribute 
to the decay, the relative phase with respect to the quarkonia components is 
at this stage irrelevant; by inverting / — > — / we would find the same results 
for the decays, but opposite glueball-quarkonia phases. In turn, the relative 
phases of the N and S components are not symmetric under e — > — e; as 
discussed above, in / (1370) they are in phase, while in / (1710) they are out 
of phase. The state |/o(1500)) behaves like a N state with a decreased width, 
while the S component is small [9,21]. Thus the decay into KK is smaller 
than the tttt channel. The iV and S state components are in phase contrary 
to the results of [5,9]. 
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The experimental uncertainties of the /o(1370) resonance are large, no average 
or fit is presented in [1]. The main problem connected with this resonance is 
its large width (200-500 MeV) and its partial overlap with the broad low-lying 
a = / (400 - 1200). However, the results from WA102 [20] indicate a large N 
component in its wave function; the results from Crystal Barrel (summarized 
in [22] and subsequently analyzed in [23]) confirm such a trend (see also [24] 
for a recent review). 

Table 2. Decays of /i = / (1370). 

Quantity Exp (WA102) Theory 

T fl ^ RK /T fl ^ nn 0.46 ±0.19 0.34 

Y h ^ m /Y fl ^ n 0.16 ±0.07 0.06 

{Y h ) 2p (MeV) "small" 166 

Predictions for the two-pseudoscalar decay modes are in acceptable agree- 
ment with the results of WA102 as shown in Table 2. The measured ratio 
Y h ^/Y fl _ n7r = 34.0ig 2 [20], although the errors are large, points to a domi- 
nant 4-7T contribution to the total width. Our prediction gives however a sizable 
contribution of the two-pseudoscalar decay mode. 



An analysis by Crystal Barrel [23] also indicates sizable partial decay widths 
of the 4tt decay channels: Y fl ^ aa = 120.5 ± 65 MeV and Y h ^ pp = 62.2 ± 
28.8 MeV. The same analysis gives the following two-pseudoscalar partial 
widths [23,24]: Y h ^ n = 21.7 ± 9.9 MeV, Y fl ^ KR = (7.9 ± 2.7 MeV) to 
(21.2 ±7.2 MeV) , Y h ^ m = 0.41 ±0.27 MeV. On the contrary, the analysis of 
[25] reports the ratio Yf 1 ^ !rn /(Yf 1 ) tot = 0.26 ± 0.09, which is in disagreement 
with the previously discussed experimental results. 

In the original work [5], a quarkonium state N has a very large two-pseudo- 
scalar width with (rjv) 2 p ~ 500 MeV. Since the /o(1370) is interpreted as 
a dominant N state [5] one also expects a large value for (Yf l ) 2 p. A sizable 
value for (rjv) 2 p is also predicted in [9]; for the mixed state / (1370) one 
has (Yf 1 ) 2 p = 115.7 MeV comparable to the present result. In a recent fit 
[7] small two-pseudoscalar partial decay widths are obtained by the following 
mechanism: for /o(1370) the glueball component gives large decay amplitudes 
which interfere destructively with the contribution of the N component. As 
a result [7], the two-pseudoscalar decay width (I7J2P is smaller than (Y f 2 ) 2P 
with {Y h ) 2P /{Y h ) 2P = 10.0±3.0; at the same time (r /3 ) 2P /(r /2 ) 2 p = 0.7±0.2. 
In the present fit with a suppressed glueball decay we have the following decay 
widths into two pseudoscalar pairs: (r/J 2 p = 165.79 MeV > (r/ 3 ) 2 p = 143.27 
MeV > (r/ 2 ) 2 p = 50.82 MeV. The experimental situation concerning the 
absolute scale of the partial decay widths into pseudoscalars remains unclear, 
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but the explicit values constitute a key test of mixing scenarios. 

The theoretical partial widths of the /o(1500) are in good agreement with the 
data (see Table 1) apart from a slight underestimate of the 2rj channel. We 
also obtain Tf 2 _ m ' = 0.036 MeV as compared to the experimental value of 
Tfc-tjrf = 2 ± 1 MeV. Taking into account the finite width of the resonance 
will lead to an increase of the theoretical value. 

For the decays of /o(1710) we summarize our results compared to the data of 
WA102 [20] in Table 3. 

Table 3. Decays of / 3 = /„(1710). 



Quantity 


Exp (WA102) 


Theory 




5.0 ±0.7 


4.70 


r '/3->W /I" '/3-»inr 


2.4 ±0.6 


1.17 




< 0.18 


1.59 


(r /3 ) 2p (MeV) 


"dominant" 


143.27 



The first two ratios, already included in the fit of Table 2, can be reproduced. 
The theoretical ratio T f 3 ^ vrj i /T f :i _^ n7T , which is not included in [1], is in com- 
plete disagreement with the WA102 result. The dominance of the rjrj' mode 
over Tin is a solid prediction, which does not depend very much on the choice 
of parameters. A confirmation of the experimental result could possibly hint 
at a sizable role of direct glueball decay. 

The ratios of two-pseudoscalar decay modes of a (1450), included in the fit of 
Table 2, are well reproduced. The prediction for the two-pseudoscalar width 
of (r„ ) 2 p = 84.26 MeV is smaller than the total width of 265 ± 13 MeV. 
However, the experimental ratio (r oo _ ( . a ,^ 7r /r ao _ > ^) is not known: no average 
or fit is listed in [1]. The experimental value from [26], which is 10.7 ± 2.3, 
would imply a dominant urni mode and in turn rather small two-pseudoscalar 
partial decay widths. This finding is in disagreement with the results of [5, 27]. 
In [5] a value of (r ao ) 2 p = 390 ± 110 MeV is found, in the work of [27] 
one has (r ao ) 2P = 420-940 MeV for a masses in the range of 1200-1400 
MeV. For our value for (r ao ) 2 p we obtain the estimate: r ao _ >a , 7r „./r oo _ >7rj; = 

[(r ao )tot - (r O0 ) 2 p]/ r O0 _^ ~ 4.5 . 

Our result for Yk*~*Ktt underestimates the experimental value by about a 
factor of 5 (see Table 1). Furthermore, for the additional Krj decay channel 
we get T K *^ Kr) /T K *^K = 0.026. 

In [27] a value of T K *^ Kw = 340 MeV is predicted, but, as discussed above, 
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(r ao ) 2P is of the order of 1 GeV, much larger than the full width. Similarly, 
in [5] with T K *^ K7T = 200±20 and (r ao ) 2P = 390±110 MeV the first result un- 
derestimates while the second overshoots the experimental value. A full anal- 
ysis in the 3 P model [28] results in V K *^ K = 166 MeV, T N ^ W7r = 271 MeV. 
Unfortunately, the resonance a (1450) is not discussed in [28]. The authors of 
[28] also tried to fix Tk*^k to its experimental value, and then calculate the 
27r partial width of a N state, obtaining Tn^ W7T ~ 450 MeV. The last result 
implies a very large two-pseudoscalar and full width for a iV state. A full ex- 
perimental determination of all relevant decay modes involving a (1450) and 
Kq(1A30) would certainly help to clarify this issue. 

If /o(1370) is dominantly nn, as in [5-9, 11], there is, as discussed above, an in- 
compatibility of the present experimental small two-pseudoscalar partial decay 
widths [23] and various model calculations. At the same time, the consistent 
understanding of the isodoublet states IT* (1450), the isovectors a (1450) is 
still incomplete. 

Proceeding as in a [28], that is fitting T k *^k (together with the ratios of 
a (1450)) we find c s d = 0.0180 GeV (larger than a factor 2 when compared to 
the previous fit in Table 1) and c s m = 0.0073 GeV. With these values one has 
T k *^k = 281 MeV, T ao ^ KR /T ao _ t7TV = 0.88 and V /V ao ^ v = 0.30 and 
^a a ^2P = 381.17 MeV. The corresponding decay width of the N state into 
two pseudoscalars is T N ->2P = 937 MeV (for a mass of M N = 1.448 GeV). 
We then find similar results as in [5,28], implying that the discussed trend 
is rather model independent. Fitting the remaining free parameters to the 
experimental data, a \ 2 minimum results in: = 1.453 GeV , Mq = 1.504 
GeV , M s = 1.698 GeV , / = 0.0281 GeV 2 , e = 0.209 GeV 2 . The mixing 
parameter / is much smaller than in the full fit, in disagreement with the 
results of [5, 7, 9]. As a further consequence the state ~ \N) has a very large 
two-pseudoscalar width of about 1 GeV. At the same time I/3) ~ \S) , which 
also has a large (and not acceptable) two-pseudoscalar width of ~ 700 MeV. 

As a further consequence we discuss the two-photon decay rates of the scalar 
resonances. We assume that the coupling cf is suppressed with respect to c s e , 
i.e. we set the glueball-photon coupling cf to zero. The ratios of radiative decay 
widths as a prediction of the fit are: 

r fl ^ 2l : r /2 ^ 27 : r /3 ^ 27 : r a o_ 27 = 1 : 0.367 : 0.004 : 0.491 , (6) 

which are independent of the coupling c s e . The result for r f 2 ^2-y/T is 
in qualitative agreement with the results of [7, 11]. The ratio r 7i->2 7 ) 
however, is considerably smaller than in the previous works. The suppression 
of originates from the destructive interference between the iV and S 

components, which in turn is traced to the flavor mixing with e > in ac- 
cord with the phases of [17,18]. Another interesting prediction is the ratio 
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r a ;|^27/T/i->27) which is relatively large. 

The experimental status of the two-photon decays is still incomplete. For the 
/o(1370) two values are indicated in PDG2000 [29]: 3.8± 1.5 keV and 5.4±2.3 
keV. However, it is not clear if the two-photon signal comes from the /o(1370) 
or from the high mass end of the broad /o(400 — 1200). The PDG currently [1] 
seems to favor this last possibility, but the data could also be valid for the 
/o(1370). We therefore interpret the two experimental values as an upper limit 
for the two-photon decay width of the / (1370). Signals for two-photon decays 
of / (1500) and / (1710) have not yet been seen; the following upper limits 
are reported [1]: 

r/ (1500)^2 7 (r/ (1500)-+7r7r/r/ (l500)tot) < 0.46 keV, 

r/ (mo)- > 2 7 (r /() (i7io)^x/r /o (i7io)tot) < O.ll keV. (7) 

Using the known branching ratio r /o(150 o)^27r/r/ (i5oo)tot one gets T fo{1500) ^ 2 y < 
1.4 keV [21]. An accepted fit for ^ f ^mo)^KK f {i7io)tot is n °t reported in 
[1]. Using the value from [30] with r /o(1710) ^^/r /o(m0 ) tot = 0.38io;?3 we find 
an upper limit of the order of Fj ( 1710 )^27 ~ 0.3 keV. 

For an absolute prediction of the two-photon decay widths we use c s e = 0.0761 
GeV -1 as determined in the model approach of Ref. [11]. For the non-strange 
quarkonium state we get Tjv^27 = 0.904 keV, while for the isovector and 
mixed scalars we have: 

Y h ^ 2l = 0.702 keV , T h ^ 2l = 0.258 keV , 

r /3 ^ 27 = 0.004 keV , r a o^ 27 = 0.345 keV . (8) 

The results for the mixed states are below the current upper limits. 

The estimate for the 27 decay of the bare quarkonium state nn of 0.904 keV 
is smaller than the one of [21], where the following expression has been used: 

r ^ (0++) = fc (^P) 3r ^ (2++) (9) 

The coefficient k is 15/4 in a non-relativistic calculation, but becomes smaller 
when considering relativistic corrections [31]. In [21] a range of values for 
k from 2 to 15/4 is considered. Our chiral Lagrangian approach combined 
with [11] points to a smaller value of k; using our result for Fjv_^ 2 7(0 ++ ) and 
taking the value T N ^ 2l (2 ++ ) = 2.60±0.24 keV [1] at M N {2 ++ ) = 1.27 GeV we 
get k ~ 0.25. This result is model dependent, since it relies on the parameters 
for the covariant description of the scalar mesons used in [11]. Hence a fully 
covariant treatment implies strong deviations from the non-relativistic limit. 
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Results when including glueball decays - In line with strong coupling argu- 
ments [5] direct glueball decay can be suppressed, leading to the phenomenol- 
ogy discussed previously. Including explicit glueball decay in a full fit to the 
data points of Table 1 with nine parameters (three bare masses, two mix- 
ing parameters / and e, four decay strengths c% c s m , c 9 d , c^J does not lead to 
a meaningful and unique x 2 minimum. To study the effect of direct glueball 
decay in the context of an effective Lagrangian, we pursue following strategy: 
we keep the bare masses listed in (4), which are in accord with other phe- 
nomenological considerations, and only fit the mixing parameters and decay 
strengths. 

We first consider the 577(3) flavor symmetry limit for the direct glueball decay, 
that is c g m = 0. A fit to the data of Table 1 results in the parameter values: 
/ = 0.0650 GeV 2 , e = 0.211 GeV 2 , c s d = 8.50 MeV, c s m = 2.58 MeV and c g d = 
—0.004 MeV with x 2 = 29.00. The flavor symmetric glueball decay constant c 9 d 
is very small in this fit, whereas the other parameters are practically unchanged 
resulting in a x 2 only sligthly smaller than before. 

Flavor symmetry breaking in the direct glueball decay is encoded in the 
strength parameter c 9 m . Inclusion of this additional parameter in a full fit 
results in the set of parameters: / = 0.0654 GeV 2 , e = 0.203 GeV 2 , c s d = 8.90 
MeV, c s m = 0.87 MeV, c 9 d = 0.42 MeV and c 9 m = -1.17 MeV with X 2 = 28.83. 
The dominant parameters /, e and c s d remain stable, when including the di- 
rect glueball decay. The coupling c s m decreases slightly, whereas the glueball 
coupling c d remains suppressed and is smaller than c g m . Latter effect might 
suggest a strong flavor symmetry breaking when discussing direct glueball 
decays. Again, the resulting masses of the scalar-isoscalar mesons and decay 
rates listed in Table 1 remain essentially unchanged for a pure fit to the data. 
Without further guidance from theory concerning the direct glueball decay the 
present phenomenological fit seems to indicate that this decay modes are sup- 
pressed in line with the original arguments given in [5]. This result, however, 
contradicts the approach chosen in [7], where direct glueball decays dominate. 

In the context of an effective chiral Lagrangian, evaluated at tree level, we 
discussed the phenomenological consequences for the scalar meson sector. 
Thereby, the scalar glueball is assumed to mix minimally with the scalar- 
isoscalar quarkonia states. A fit of the mixing parameters and coupling streng- 
ths entering in this scheme to the accepted data yields several results: the 
direct flavor mixing between quarkonia states [17] leads to a structure of the 
physical f states, which differs from the conventional scheme as in Refs. [7, 9]. 
The different pictures can be directly tested in radiative 27 decays of the / 
mesons. The relevance of direct glueball decays, which play a dominant role 
in [7], cannot be verified in the present scheme. Further input from lattice 
simulations, as for example demonstrated in Ref . [32] , could serve to obtain a 
better understanding of this issue. 
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